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Abstract 

We compute in the instanton vacuum the nucleon matrix elements of the twist-4 
qv . QCD operators describing power corrections to the second moments of the unpolar- 

ized structure functions, Fl and F2. Our approach takes into account the leading 
contribution in the packing fraction of the instanton medium, p/R <C 1. Paramet- 
rically leading are the matrix elements of a twist-4 quark-gluon operator, which 
are of the order of the inverse instanton size, 1/p 2 = (600 MeV) 2 . The matrix el- 
ements of the four-fermion (diquark) operators are suppressed by a factor (p/R) 4 
and numerically small. These results are in agreement with the pattern of phe- 
nomenological l/Q 2 -corrections to R = ctl/&t and F% found in QCD fits to the 
data. In particular, the rise of it! at low Q 2 can be obtained from instanton-type 
vacuum fluctuations at a low scale. 
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1 Introduction 

An outstanding problem in the theory of deep-inelastic scattering is to understand the 
transition from the asymptotic region, where the scale dependence of structure functions 
is described by perturbative evolution, to the region of low Q 2 , where non-perturbative 
effects play an essential role. Starting from high Q 2 , the onset of non-perturbative effects 
manifests itself in power (1/Q 2 -) corrections, which are governed by non-perturbative 
scale parameters. Aside from the well-understood target mass corrections there are dy- 
namical power corrections, which in QCD are determined by nucleon matrix elements of 
operators of non-leading twist (> 2) |j], @,@|. In the partonic language, they describe the 
effects of the interaction of the quarks with the non-perturbative gluon field in the nu- 
cleon, or of 2-body correlations between quarks, on the DIS cross section. A well-known 
example are the longitudinal structure function F^ = F 2 — 2xFi, or the ratio R = Cl/c^t, 
which both happen to be zero in the naive parton model (i.e., without radiative correc- 
tions). QCD fits of the data, using existing NLO parametrizations to describe the twist-2 
contribution, leave room for sizable l/Q 2 -power corrections [|], j§, |§, 0. It remains a 
challenge to understand the origin of these corrections from first principles. 

During the last years there has been mounting evidence for an important role of 
instanton-type vacuum fluctuations in determining the properties of the low-energy ha- 
dronic world ||. The so-called instanton vacuum is a variational approximation to the 
ground state of Yang-Mills theory in terms of a "medium" of instantons and antiinstan- 
tons, with quantum fluctuations about them, which is stabilized by instanton interactions 
|J. The strong coupling constant is fixed at a scale of the order of the inverse average 
size of the instantons in the medium, p _1 ~ 600 MeV. The most important property of 
this picture, which had been postulated on phenomenological grounds by Shuryak [|Kj 



and was later established in the variational calculation of Diakonov and Petrov || , is the 
diluteness of the instanton medium, i.e., the fact that the ratio of the average size of in- 
stantons in the medium to their average distance, R, is small: p/R ~ 1/3. The existence 
of this small parameter makes possible a systematic analysis of non-perturbative effects 
in this approach. 

In particular, the instanton vacuum explains the dynamical breaking of chiral sym- 
metry, which is crucial in determining the structure of strong interactions at low energies 



TT| , |12f . It happens due to the derealization of the fermionic zero modes associated with 
the individual (anti-) instantons in the medium. Using the l/iV c -expansion one derives 
from the instanton vacuum an effective low-energy theory, whose degrees of freedom are 
pions (Goldstone bosons) and massive "constituent" quarks. The nucleon is described 
as a chiral soliton in this effective theory, as a state of N c "valence" quarks bound by a 



classical pion field JIB 



The instanton vacuum has also received strong support from direct lattice simulations. 
Not only were the basic characteristics of the instanton medium reproduced after "cooling" 
of the quantum fluctuations ||14]| , also the correlation functions of various hadronic currents 
were found to remain practically unchanged under cooling, showing that the instanton 
vacuum takes into account the relevant non-perturbative fluctuations in these channels. 



Given the general success of the instanton vacuum in describing hadronic properties, it 
is natural to ask if this scheme of approximations can explain the non-perturbative effects 
observed in the scale dependence of structure functions at low Q 2 . More ambitiously, 
one could also ask the opposite question: Can we use the rich and increasingly accurate 
information obtained from DIS to learn something about properties of non-perturbative 
vacuum fluctuations? 

The aim of this paper is to study the l/Q 2 -power corrections to unpolarized struc- 
ture functions, in particular, to the longitudinal structure function, F L , in the instanton 
vacuum. We restrict ourselves to the lowest moments, for which the power corrections 
can be studied using the operator product expansion (OPE) |Ij, @- In particular, we use 
the instanton vacuum to compute the nucleon matrix elements of the twist-4 operators, 
both the quark-gluon and the four-quark operators, and discuss their relative magni- 
tude. Our approach takes into account the leading contribution in the packing fraction 
of the instanton medium, p/R |T5|, [Tj]j . Our main conclusion is that the matrix element 



of the twist-4 quark-gluon correlator is parametrically large and accounts for most of the 
observed scaling violations in Fl, while the four-fermionic operators are suppressed. 

The quantitative description of power corrections is a complex problem. The set 
of higher-twist operators used in the OPE is over complete, since matrix elements of 
different operators are related by the QCD equations of motion. In order to have results 
which are independent of the choice of the operator basis one needs to make consistent 
approximations in dealing with the QCD operators and with the nucleon state, something 
which is generally impossible to achieve in phenomenological approaches such as the bag 
model H]. In this respect the instanton vacuum with its inherent small parameter, p/R, 



offers some decisive methodological advantages. As was shown in Ref.||16|, in this approach 
the QCD equations of motion are preserved at the level of effective operators, which is 
a consequence of the one-instanton approximation justified in leading order of p/R, and 
of the fact that the fermions couple to the instanton through the zero mode, see [0 for 
details. It is crucial in this context that the description of the nucleon as a chiral soliton 
is a fully field-theoretic one and does not involve any additional dynamical input besides 
the chiral symmetry breaking induced by the instanton fluctuations. Also, it should be 
pointed out that, contrary to other field-theoretic approaches such as QCD sum rules or 
lattice calculations our approach is based on a small parameter, p/R, which endows the 
a priori structure-less set of higher-twist matrix elements with a parametric order, thus 
greatly simplifying the discussion of power corrections. 

The plan of this paper is as follows. In Section |2|, we summarize the QCD predictions 
for the scale dependence of the unpolarized structure functions Fl and F 2 to order l/Q 2 , 
including twist-2, target mass and dynamical twist-4 corrections. The central part of the 
paper, Section [| is devoted to the calculation of the relevant twist-4 matrix elements 
in the instanton vacuum. Subsection EO describes the basic properties of the instanton 



vacuum and the resulting low-energy effective chiral theory. In Subsection |3.2| we outline 
the method by which the QCD twist-4 operators can be mapped onto effective operators, 
whose matrix elements can be computed in the low-energy effective theory [TJ|. As a 
novel feature, we show here that the QCD twist-4 quark-gluon operators which produce 



parametrically large matrix elements in p/R can be represented in the effective theory by 
one-body operators built from the massive "constituent" quark field. This representation 
turns out to be very convenient both for discussing general properties of instanton-induced 
higher-twist effects, as well as for computing the nucleon matrix elements. In Subsection 
|3.3| we compute the nucleon matrix elements of the effective operators, within the large- 
N c picture of the nucleon as a chiral soliton. In Section ^ we use the instanton results for 
the twist -4 matrix elements to estimate the scale dependence of F^ and F 2 , and compare 
with the data for R = ctl/ct and F 2 . We also compare our results for the matrix elements 
with those of other theoretical approaches. Conclusions and an outlook are presented in 
Section ||. 

Appendix |A] contains a short derivation of the QCD expressions for the dynamical 
twist -4 corrections to unpolarized structure functions using the operator product expan- 
sion m . Since calculations of higher-twist effects are highly dependent on conventions we 
explicitly state all relevant conventions there. 

We are not considering here power corrections to DIS structure functions resulting 
from instanton contributions to the coefficient functions of the OPE, i.e., semiclassical 
corrections associated with single instantons of size p ~ 1/Q |I7|, |18| . These are suppressed 



by very large powers of 1/Q and typically negligibly small. Rather, we use a medium of 
instantons of average size p = (600 MeV) -1 as a model for the non-perturbative vacuum 
structure at the hadronic scale, in order to calculate the matrix elements parametrizing 

1 /Q 2 -corrections. 

The power corrections to Ft and R = <Jl/&t have been estimated previously in a vari- 
ety of models. In the so-called transverse basis of twist-4 operators the l/Q 2 -corrections 
to Fi can be related to an "intrinsic" transverse momentum of the partons 0. This 
language is frequently used to model dynamical twist-4 contributions. A model interpo- 
lating between large Q 2 and the limit Q 2 — > has been proposed in Ref.|19|. In Ref. [pOfl 
the leading-order twist-4 contribution at large x has been related to the derivative of 
the twist-2 distribution. Also, the magnitude of the twist-4 contribution to Fl has been 
estimated from the renormalon ambiguity of the perturbation series for the 0(as) twist- 

2 coefficient functions calculated in the large- Nf limit |2TL E2|. In Ref.[|23| an attempt 
was made to model the non-perturbative twist-4 contribution at large x by introducing 
a coupling to vacuum condensates. Higher-twist contributions to unpolarized structure 
functions have also been studied in the bag model p4[ ^5[ . 



2 Power corrections to unpolarized DIS 

The hadronic tensor for unpolarized electromagnetic scattering off the nucleon has two 
independent symmetric structures, which are parameterized in terms of the two structure 
functions F 1 and F 2 , or, equivalently, F L = F 2 — 2xF x and F 2 : 



W, v { V ,q) = ^Jd'ze^iplJ^JMlp) 



F L (x,Q 2 ) ( _ g^ 
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+ F 2 (x,Q 2 ) 
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Here, J^O 5 ) = ^K- 2 ) ^ s the electromagnetic current operator, p^ the nucleon momentum 
(p 2 = M^), x = —q 2 /{2p-q) and Q 2 = — g 2 . Here and in the following, all matrix elements 
are understood to be averaged over the target spin, 

(p\...\p) ee lJ2( P S\...\pS), (2.2) 

and the states are normalized according to (p'\p) = 2p°(2ii) 3 S^(p' — p). In connection 
with experiments one frequently uses a slightly modified definition of the longitudinal 
structure function ||: 

F L (x,Q 2 ) = F 2 (x,Q 2 )(l + ^^-)-2xF 1 (x,Q 2 ) 



O 



A M 2 T 2 

x,Q 2 ) + ^-F 2 (x,Q 2 ). (2.3) 

This function is directly related to the experimentally measured ratio R = ai/a?, namely 

In QCD, the electromagnetic current is carried by the quark field, 

W = Y, e rff(*) , rM*)> ( 2 - 5 ) 

where e/ are the quark charges. The scale dependence of the moments of the structure 
functions up to level 1/Q 2 can be determined by relating the hadronic tensor to the imag- 
inary part of the forward amplitude for the scattering of a virtual photon off the nucleon, 
and applying the operator product expansion (OPE) to the time-ordered product of elec- 
tromagnetic currents in the latter |l|, H|. A short summary of this approach, including all 
relevant conventions, is given in Appendix |A|. We consider here the second moments of 
the structure functions: 

^■<M = r dx { p ^ . (2 ,) 

M 2 (2,Q 2 ) J A \ F 2 (x,Q 2 ) 

We distinguish between the scaling (non-power suppressed) part of these moments, which 
comes from the contribution of twist-2 spin-2 operators in the OPE, the target mass 
corrections (oc M^/Q 2 ), which arise from twist-2 spin-4 operators, and the dynamical 
power corrections coming from contributions of operators of twist 4 and spin 2. 



In the case of the longitudinal structure function the twist-2 contribution is of order 
as, since Fp is zero in the naive parton model due to the Callan-Gross relation. The 
coefficient functions have been computed in Ref.[||, see also Refs. [^, |2B|. In leading 
orderQ 
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where Cp = 4/3 and the sums run over all light flavors (u, d and s) in the target. Here 

r( n ) r(«) ■ 
^f ) -^g c 

the target, 



Lj , Lq denote the matrix elements of the twist-2 spin-n quark and gluon operators in 



(p\ ^/7 M i^ V M2 • • • W^ipf \p) - traces = 2L ( , n) 



J f 



x (P w ••• JW - traces) , (2.9) 
Mn i^/ - — <g 

where on the L.H.S. complete symmetrization in \i\ . . . fi n is implied. Here 



(p\F, lQ D, 2 ...D, n _ 1 F^ n \p) -traces = 2L^ 



d,-iA;(x)^, 



(2.10) 



with tr[A a A 6 ] = 2S ab , is the covariant derivative in the fundamental representation of 
SU(3), 



flV 



* [V M , V,] 



F a — 
nv 2 



the gauge field, and 



(n F Y = 8 F a - f abc A b F c 



(2.11) 



(2.12) 



the covariant derivative in the adjoint representation. 

The target mass corrections to the second moments are proportional to the matrix 
element of the spin-4 twist-2 operator (see Refs.|I|, |3[ and Appendix ^): 



M L (2,Q 2 ) targetmass 
M 2 (2,Q 2 ) targetmass 



Q 2 



/ 
/ 



2 Q2^ e f L f 



(2.13) 
(2.14) 



■"■We write down the expressions for Ft,, Eq.(2.3); the only difference to the original Fl, will be in the 
target mass corrections. 



where we have restricted ourselves to the 0(a.g) (tree-level) coefficient function^. 

Our main object of interest are the dynamical power corrections, which are determined 
by matrix elements of twist-4 spin-2 operators |TJ, |2|]. There are two types of twist— 4 
contributions, shown schematically in Fig.[l], resulting from different contractions of the 
quark fields in the two electromagnetic current operators. Contribution (a) describes the 
interference of the scattering of the photon off a free quark with scattering off a quark 
interacting with the non-perturbative gluon field in the target; this contribution gives rise 
to matrix elements of twist-4 quark-gluon operators in the target. Contribution (b) de- 
scribes the interference of scattering off two different quarks in the target, accompanied by 
a perturbative gluon exchange; it leads to matrix elements of four-fermionic ("diquark") 
operators in the target. In the longitudinal structure function the latter contribution is 
absent. We have recalculated the coefficient functions of the twist-4 operators, following 
the approach of Ref . [[J , and reproduced the results given in that paper. A short outline 
of the calculation is given in Appendix |A]. For the second moment of the longitudinal 
structure function one finds 

M L (2,Q 2 ) twist - 4 = ^E4(-^/ + i 5 /)- ( 2 - 15 ) 

Again we have restricted ourselves to the tree-level coefficient function. Here, Af and 
Bf [their dimension is (mass) 2 ] denote the reduced matrix elements of the twist-4 spin-2 
operators: 

A a /3,f = ^/7 Q! (-D 7 F 7/ g)V'/ - trace, (2.16) 

B a p,j = ^/7 7 75« (Vo.F/37 + F ai V^ i/jf - trace, (2.17) 



with 



A a/3,f I . . n [ Af 



(p\i B >\p) = 2< B \ x {p^-\p 2 g a p). (2.18) 

It is understood that the operators are symmetrized in the Lorentz indices according to 
A a /3 — ► \(A a p + Ap a ); this will not be explicitly indicated in the following. Here 

W v = -e^ pa F pa (2.19) 

is the dual field strength. We are using the convention 75 = — Z7°7 1 7 2 7 3 (as in Refs.JT], |ll], 
IR |TB[), which differs from the one of Bjorken and Drell by a minus sign, and e 0123 = 1. 



Note that, using the QCD equations of motion, one can replace the forward matrix element 
of the operator A a @^ f by that of a four-fermionic operator: 

Aa ,f - -g 2 Uf\i^A (E#y7/3#j- (2.20) 



For the second moment of Fl instead of i7 the coefficient in Eq.(2.13) would be —3 instead of +1, 
cf. Eq.Q. 



In the case of the structure function F<i one has to consider also the twist-4 contri- 
butions of type Fig.|l] (b), which can be computed using standard Feynman rules. One 
finds 

l^o/5. 1 \ 1 



M2( 2 5g2)tw ist-4 = —Y,e}(--A f --B f ) + —Y,e f e r C ff , 1 (2.21) 



where Cff denotes the reduced matrix element [cf. Eq. (|2.18 )] of the four-fermionic op- 
erator 

^/y7 a 75^/J (^/'y7/375^/'j • (2.22) 

3 Twist— 4 contribution to Fi and F<i from instantons 

3.1 Instanton vacuum and low— energy effective theory 

The framework for our calculation of matrix elements of QCD operators is the instanton- 
based vacuum of Diakonov and Petrov ||, which is a variational estimate of the Euclidean 
Yang-Mills ( "quenched QCD" ) partition function in terms of configurations consisting of 
a superposition of instantons and antiinstantons (J's and J's), with quantum fluctuations 
about them. The medium of J's and J's stabilizes due to an effective repulsion between 
the pseudoparticles. The resulting statistical ensemble turns out to be well described by 
a partition function of independent J's and J's with an effective size distribution, with 
average size p = (600 MeV) -1 . In the large-iV c limit the size distribution is sharply 
peaked, so one may take the sizes of all J's and J's to be equal to p. An important 
property, which is in fact crucial for this picture to be consistent, is the diluteness of the 
medium of J's and J's. In the approach of Ref. || the ratio of the average size of instantons 
to their average distance in the medium was found to be p/R ~ 1/3 <C 1. This fact can 
be traced back to the "accidentally" large value of the coefficient of the one-loop beta 
function of QCD, b = 11NJ3 « 10. 

When quarks are included, the dynamical breaking of chiral symmetry is a conse- 
quence of the interaction of the quarks with the fermionic zero modes associated with the 
individual instantons [11 1. The most compact way to express this is through an effective 



low-energy theory, which has been derived in the large- N c limit using the so-called zero 



mode approximation (which is exact in leading order of p/R) [|12| |15| . It can be formu- 
lated in terms of quark fields (Nf flavors) with a dynamical mass and a Goldstone bosons 
(pion) field, with an effective action 

N f <- -> 

S eff = fd 4 x ]T $,{x) <V/HbA "- tMF(-d 2 )U^(x) rf F(-d 2 ) i/> f (x). (3.1) 

/,/'=! L J 

We have passed to Euclidean field theory by continuing the time dependence of the fields 
to imaginary times and introducing Euclidean fields and gamma matrices according to the 



conventions given in Appendix |B|. (Unless stated otherwise, all formulas in Subsections [TT] 
and |3.2| will be written in terms of Euclidean quantities.) In Eq. (|3.1| ), M is the dynamical 
quark mass; parametrically it is of order 



Mp 



R 



The pion field, 7r(x), couples to the quarks in a chirally invariant way 



U^(x) 



±¥*u(z) + l - 



x 



U(x) 



exp [iir a (x)T a ] 



(3.2) 



(3.3) 



Furthermore, F(—d 2 ) is a form factor proportional to the wave function of the fermionic 
zero mode of the instanton in momentum representation [|TT|, [HJ , 



F(k 2 



d_ 
dt 



'- ' I) K " G 



"(l)**G 



pVP, (3.4) 



where I n and if n are modified Bessel functions. It satisfies F(0) = 1, and F(k 2 ) ~ 
(p 2 k 2 )~ 3 / 2 for fc 2 — > oo, i.e., F{k 2 ) drops to zero for momenta of the order of the inverse 
instanton size. Here and in the following we use a shorthand notation for the action of 
the form factors, 



F(-d 2 )ip{x) 
^(x)F(-d 2 ) 



J^k_ e±lHx -, )ne) | 



ij(x') 



(3.5) 



Thus, the inverse instanton size plays the role of an ultraviolet cutoff of the effective 
theory. It is important to note that the way in which the ultraviolet regularization of 
the effective low-energy theory is implemented, namely through form factors restricting 
the virtuality (Euclidean momentum) of the quark fields, is unambiguously determined 
by the instanton vacuum. 

The effective action, Eq. fl3.lp , has been derived in the large-iV c limit, and can be used 
to compute hadronic correlation functions in a 1 /^-expansion (saddle point approxi- 
mation) (12| . Correlation functions of mesonic currents computed in this approximation 
generally show good agreement with phenomenology [TT, [T2J. The nucleon correlation 
function at large iV c is described by a saddle point with a non-trivial classical pion field 
( "soliton" ) [I3| . In the nucleon rest frame it is static and of "hedgehog" form, 



^class(x) 



exp 



-P(r) 



r = 



(3.6) 



where P(r) is called the profile function, with P(r) — > for r — > oo. Integration over 
translational and (iso-) rotational zero modes of the saddle point field gives rise to nucleon 
states of definite spin/isospin quantum numbers [H|, with the delta resonance appearing 
as a rotational excitation. This picture of the nucleon gives a good description of practi- 
cally all hadronic observables of the octet and decuplet baryons (for a review see Ref . 1 27| ) . 
It also describes well the twist-2 parton distributions of the nucleon 



10 



3.2 Effective operators for QCD twist— 4 operators 

The instanton vacuum is a variational approximation to the full QCD partition function 
and thus allows to evaluate hadronic matrix elements of QCD operators, including such 
involving the gluon field. It is understood that the QCD operators are normalized at 
the scale of the inverse instanton size, p ~ p _1 , the scale at which the strong coupling 
constant is fixed when determining the properties of the instanton vacuum. 

A convenient method for computing hadronic matrix elements of QCD operators, in 
the same spirit as the effective low-energy theory, Eq.( |3.1| ), is the method of effective 



operators |T5) . Using the same approximations as those which went into the derivation of 
the low-energy theory — the diluteness of the instanton medium, the zero mode approx- 
imation, and the large-iV c limit — it is possible to "translate" QCD composite operators 
into effective operators, whose hadronic matrix elements can be computed within the 
effective theory. This method has been used to calculate various vacuum condensates 



30(1 as well as spin-dependent and -independent nucleon matrix elements of higher-twist 



operators [|H|. It was shown in Ref.[16|] that basic relations between matrix elements 
of higher-twist operators following from the QCD equations of motion are preserved at 
the level of effective operators. We now compute in this approach the twist-4 operators 
arising in the l/Q 2 -corrections to F L and F 2 , Eqs. (|2T6|) , (|2T?|) and (ggg) . 



For a QCD composite operator made purely from quark fields the effective operator 
in leading order of the packing fraction, p/R, is simply given by the QCD operator 
with the fields replaced by the massive quark field of the effective low-energy theory 
fl~5|j . For a gluon or quark-gluon operator, one must in addition (after passing to the 
Euclidean theory) replace the gluon fields in the operator by the field of one I (I) and 
integrate over its collective coordinates; multi-instanton contributions are suppressed by 
additional powers of the packing fraction. One immediately sees that at this level the 
effective operator for the QCD operator A a p, Eq.( [2.16| ), is identically zero, since the 
instanton field is a solution to the Euclidean Yang-Mills equations, 

D 7 F 7a (x) J(i) = 0. (3.7) 

Thus, the operator A a p requires at least two-instanton contributions to be non-zero, and 
its matrix element is of order (p/R) 4 , i.e., parametrically suppressed. In contrast, the 
function of the gauge field appearing in the operator B af3 , Eq.( ^.17[ ), is non-zero at one- 
instanton level (see below), and its matrix elements can be of order unity in the packing 
fraction. 

One may ask what would have happened if we had started instead from the four- 
fermionic version of the QCD operator, A a p, Eq.( [2.20| ), which is equivalent to the original 
quark-gluon operator by the QCD equations of motion. In Section FJ71 we shall explicitly 



compute the nucleon matrix element of the four-fermionic operator Eq. ( [2.20|) and see that 



it, too, is suppressed by a power of the instanton packing fraction: the matrix element 
is of order M 2 , while that of the operator B a p is of order pr 2 . Either way, we arrive at 
the same conclusion: In the instanton vacuum the quark-gluon operator proportional to 
the covariant divergence of the gauge field, A af3 , Eq. ([2.16|) , is suppressed relative to the 
operator B a p. 

11 



Also in Section [T4| we shall investigate the matrix elements of the four-fermionic 
operator, C a/ 3, Eq. (|2.22|) , which appears in the power corrections to F 2 only. Its matrix 
elements also turn out to be suppressed relative to those of the quark-gluon operator, 

We now compute the effective operator for B a p, Eq.( ^.17| ), in the instanton vacuum, 

For use with the Euclidean theory we 



following the steps described in Refs.]!^, [R 

define Euclidean vector components of the Minkowskian operator, Eq.( |2.17| ), according 

to (i,j = 1...3) 



(£, 



(B, 



ij)Mi 



(B 



44 E 



-(B 



00 ) M, 



(Bi4>l 



[Bio) M , (3.8) 



and express them in terms of the Euclidean fields and gamma matrices, using the con- 
ventions given in Appendix [B]. (We suppress the labels denoting Euclidean components 
in the following; unless stated otherwise all formulas in this section are for Euclidean 
objects.) It will be convenient in the following to separate the parts originating from the 
ordinary derivative and the gauge potential in the covariant derivative in Eq.( ^.17[ ): 



B, 



<*P,f 



X 



X" 



-#}(a;)— Jj^ida ^ f (x) 



<fr(.X) 



X b X c 



+ H)V,}(*h 7 75 \ T , T \ V7(s) A b a (x)F^(x) 



2 ' 2 



(3.9) 
(3.10) 



(symmetrization and trace subtraction will not be explicitly shown in the following). 
This separation refers explicitly to the "singular" gauge of the instanton field. In the first 

part, Eq. (p.9|), we have performed an integration by parts (d=d — d), dropping total 
derivatives which do not contribute to the forward matrix element. Following Refs. [p~5| , |16f 
we have to substitute in Eqs. (|3.9| ) and (|3.10|) the gauge field of an I (I) and integrate over 
its collective coordinates. For the part Eq. ( |3.9|) we need the dual field strength, which for 
/ (/) in standard orientation and centered at zero takes the form 



F a M {x) m 

JlJ,v,p~f\X) 



± *&(*)/(/) 



±(v T ) a uuf^M x )' 



iff 



(x 2 + p- 



>)2 1, 



/ X^X/3 



ar 



S-, v + 



Jr yv 



X 2 



li/3 



-iPiipv^v 1 



(3.11) 



where [rj^Y = f[^ v ,r(^ v are the 't Hooft symbols. For operators of this type the form 
of the effective operator has been derived in Ref.Jl6|. After integration over instanton 
coordinates one obtains from the part Eq. (|3.9|) the following Euclidean effective operator: 



a a/3,f 

%M 



x 



X a 



i^}{x)—7j75id a ipf(x) 



x 



iM r Nf 

— J d z J2fwPy( x - z ) 

c 9,9' 



A" 



*) F{-d 2 ) -^75 V£{z) F{-d 2 ) i, gl {z) 



(3.12) 
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Here a^ = (i/2) [7^, 7„] . The form of the effective operator is intuitively plausible: The 
gluon field in the QCD operator, Eq. (|3.9|) , has been replaced by the field of a single 
/(/), which couples to the fermions through the zero modes. This induced fermion vertex 
is chirally odd (oc a^ v ) and accompanied by a factor M; however, chiral invariance is 
maintained due to the presence of the pion field. 3 The presence of 75 in the instanton- 
induced vertex in Eq. (|3.12| ) is a consequence of adding instanton [oc (1 + 75) /2] and 
antiinstanton [oc (1 — 75) /2] contribution with different sign, cf. Eq. Q3.llp . The effective 
operator can graphically be represented as in Fig|J (b). 

In a similar way one can derive the contribution of the part Eq.( [3.10|) to the effective 



operator. The anticommutator of Gell-Mann matrices contains a color-singlet as well as 
a color-octet piece. We shall not write down these contributions here. It turns out that 
the "derivative" part, Eq. Q3.9p , gives the (by far) numerically dominant contribution to 
the matrix element, so we shall drop the contributions from Eq. (|3.10|) . 

The effective operator, Eq.( |3.12| ), can now be inserted in correlation functions of 
hadronic currents computed within the effective low-energy theory. We are interested 
only in the leading contributions in p/R, viz. Mp to the hadronic matrix element. For 
dimensional reasons, the instanton-induced vertex in Eg. ( p. 12 ) is proportional to a pos- 



itive power of the instanton size, p 2 . A parametrically large contribution to correlation 
functions can thus only come from diagrams in the effective theory containing "quadrat- 
ically divergent" loop integrals, i.e., integrals giving rise to a factor of p~ 2 , which cancels 
the factor p 2 incurred from the instanton-induced vertex (remember that p also plays the 
role of ultraviolet cutoff of the effective theory) . One can easily see that such quadratic 
divergences can only arise from diagrams of the type shown in FigfJ (c), where the quark 
propagator in the loop can be taken as the free quark propagator, (— ry • d)" 1 ; taking into 
account the coupling of the quarks in the loop to the pion field would lead to contributions 
of higher order in Mp. Thus we may replace the four-fermionic operator Eq. (|3.12| ) by a 
one-body operator representing the sum of the two contractions of Fig.||| (c). Due to the 
"separable" form of the instanton-induced vertex it is of the form: 



lM - -.t^ r«, 



2 



/' 



£'V/0r)one-bod y = —J2^ f (x)T^(id)U]},(x)F(-d 2 )^ f (x) 



+ ^),{x)F{-d 2 )U]f f {x) T%(-id) V7(ar) 



. (3.13) 



The vertices T a l and T a l can be determined by computing the loops in the diagrams of 
Fig.|2] (c) in momentum representation: 

r !>) 1 , r d*k 



r (2) M [ ^)/(^)iW*0 



The vertex involving the pion field in Eq.(|3.12) is obtained from the original Nf x ./Vy—fermionic 



instanton-induced vertex by making use of the saddle-point condition of the effective theory at large N c ; 
see iFiql for details. 
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where 



F\(p - kf 

(p - w 

F[(p + k) 2 
(p + k) 



-1p + k) a (p - k) 6 [777575^^75] 



i-\-2p -k) a (p + k) 6 [0-^75757775] 



(3.14) 



Jij,u,l3^{k) 



(X X 6 J ixvfi^ \X ) 






?G{k 2 



k 



2 u ~i v 



k 2 



M/3 



\S^8 1V , 



(3.15) 



Q{k 2 



32tt' 



MW>+(MW>-? 



t = pVfc 2 . (3.16) 



Since the instanton-induced vertex in Eq.( ^.12| ) is oc a^, the vertex in the one-body 
operator, Eq.( |3.13|) , is also chirally odd and can be reduced to Dirac structures 1 and 
a K \. Note that, again, chiral invariance is preserved by the presence of the pion field in 
Eq. (|3.13|) . Computing the loop integrals Eq.( |3.14| ) we find that both structures occur 
with the same form factors. The resulting one-body operator is given by 



t B U afi,t(x] 



one— body 



M "' 



2 fy(x) {d a d p -1 d a d K a KP ) B(-d 2 )Uf f ,(x)F(-d 2 ) V/ 



' x 



+ ij)),{x) F(-d 2 ) U]f f {x) B(-d 2 ) (dj p +1 dJ K a KP ) ij} f ix) 
where B(p 2 ) is a dimensionless scalar form factor defined as 



B{p 2 



d A k 
19 k 



Q{k 2 



F[(p-kf 
(p-k) 2 



Ik 2 



Ik 



~-Cf (cos*) + U^2 + 1 <#W) - ^-Cf (cos< 
lip \2 p z I Sp 



(3.17) 



(3.18) 



Here the Cf^> denote the 4-dimensional spherical harmonics (Gegenbauer polynomial of 
index 1): 



C£ 1J (COB0) 



sin(n + 1)# 
sin# 



cos# 



k ■ p 



(3.19) 



The integral in Eq. (|3.18|) is of order p 2 ( "quadratically divergent"), as a result of which 
B(p) is of order unity in Mp ~ (p/R) 2 . Numerical evaluation of the integral gives 



B{p 2 = 0) = 0.7£ 



(3.20) 
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and shows that B(p 2 ) ~ (p 2 p 2 )~ 5 ^ 2 for large p 2 . 

To summarize, the parametrically large contributions to the matrix element of the 
quark-gluon operator B a/3 , Eq.( |2.17| ), are contained in the "quadratically divergent" part 
of the contractions Fig|J (c), which can be represented by a chirally odd one-body quark 
operator in the background of the pion field. 

The effective one-body operator, Eq.( p.l7|) , can also be written in a different form. 



Making use of the equations of motion of the effective low-energy theory, cf. Eq.( |3.1| ), 

N f <- -> _ 

- iM ]T 4 ■ ■ ■ F(~ d 2 )U]},F(- d 2 ) fy, = $ ... * 7 „ d, fy, 

f 

N f <- -* _ 

-iM $>}, F{- d 2 )U]f f F(- d 2 ) . . . i/,, = ^} (-07k d, ...i/> f (3.21) 

/' 

[the ellipses denote additional functions of d or d, respectively], we can eliminate the 
pion field in Eq.( [3.17| ) and write the operator in a manifestly chirally invariant form. 



Integrating by parts and simplifying the product of gamma matrices, we finally obtain an 
operator 

tt B» a0t f = ^\F- 1 {-d 2 )B{-d 2 )d 2 id {alp} i }f -- trace. (3.22) 

The meaning of the "inverse" form factor here is clear in momentum representation, cf. 
Eq.( |3.5|) . In this representation of the effective operator the effects of the gluon field 
of the instanton are contained in additional contracted derivatives acting on the quark 
fields relative to the spin-2 twist-2 operator, Eq. (|2.9| ). The operator Eq. (|3.22| ) measures 
the "average virtuality" of the quark in the nucleon. This is what one should expect on 
general grounds in QCD: By the QCD equations of motion, the explicit gauge fields in 
twist-4 operators in the "collinear basis" can be traded for transverse derivatives of the 
quark fields, which in turn are related to the virtuality of the quark in the nucleon |1 . 

For momenta of the quark fields p 2 <^L JT 2 we can neglect the momentum dependence 
of the form factors in Eq. (|3.22[) , so that the effective operator reduces to a local operator: 

u B n a0,f « B(0) i)\d 2 id {a ~ip } i) f - trace. (3.23) 

Upon returning to the Minkowskian theory (see Appendix |B|) this operator becomes 

"£" a/3, /[Minkowski = B(0) $f 2 id {a jp } ipf - trace. (3.24) 

It is interesting to note that the reduced matrix element of this operator in a "constituent" 
quark state [cf. Eq. fl2.18p l, *- e -> an on-shell quark state with (^Minkowski = M 2 , has a 



negative value (for any flavor /): 

5 quark = -M 2 B(0) <0. (3.25) 

However, as we shall see below, in the nucleon the flavor-singlet matrix element is posi- 
tive, since the quarks in the nucleon have on average spacelike momenta, so (^Minkowski 
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effectively becomes negative. This exercise shows that the effective operator, Eq.( |3.22| ). 
should not be used out of context: It refers explicitly to the effective low-energy theory 
derived from the instanton vacuum, where the quarks inside hadrons have virtualities 
(Euclidean momenta) up to the ultraviolet cutoff, p~ x . In particular, the notion of "av- 
erage virtuality" of the quarks has no absolute meaning but is defined by the ultraviolet 
regularization of the effective theory. 

We remark that for the contribution from the part Eq.(RTTDI) to the effective operator 



one can also derive an effective one-body operator similar to Eq.( |3.17| ). Explicit calcu- 
lation shows that for this part the form factors are numerically an order of magnitude 
smaller than in Eq.( |3.17| ), so we can drop this contribution. 

3.3 Nucleon matrix element of the quark— gluon operator B 

We now compute the nucleon matrix element of the operator B a ^, Eq.( |2.18D , in the 
effective chiral theory. To obtain an explicit expression for the reduced matrix element, 
Bf, we contract both sides of Eq. fl2.18|) with a Minkowskian light-like vector, n^ (n 2 = 
0, rc% = M N ): 

B f = ^pn a n" <p\BafiM- ( 3 - 26 ) 

This expression is manifestly covariant and can therefore be evaluated in any frame. 
In our large-iV c approach it is convenient to work in the nucleon rest frame, where the 
nucleon is described by a static classical pion of the form Eq. fl3.6p , subject only to adiabatic 
rotations and translations whose quantization gives rise to the nucleon quantum numbers. 
We are interested in the parametrically leading contribution to the matrix element in 
the instanton packing fraction, p/R, which can be computed passing to the Euclidean 
theory and replacing the operator B a p by the effective operator obtained from instantons, 
Eq. (lCTD . 

A well-developed technique exists for computing matrix elements of quark operators 
within the l/iV c -expansion. It can directly be applied to matrix elements of the instanton- 
induced effective operators derived from twist-4 QCD operators, see Ref. |16[ for details. 



For simplicity we consider the effective low-energy theory with two light quark flavors, 
Nf = 2. The isosinglet part of the matrix element Eq. fl3.26|) appears in leading order of 
the l/7V c -expansion, the isovector part only in next-to-leading order, i.e., after expanding 
to first order in the angular velocity of the soliton. Noting that the form factor B in the 
one-body operator is of order unity in N c , standard N c counting tells us that 



B s ee Bl + Bl = 


K + B n d ~ 1, 


1 NS ee Bl-Bl = 


-{Bl-B« d ) ~ 1/N t 



(3.27) 

(the isoscalar matrix element is of the same order in N c as the momentum fraction of 
quarks, which is order unity). Incidentally, this means that the target mass corrections, 
Eqs.( p.l3] ) and (|2.14j) , are of the same order in 1/N C as the dynamical twist-4 contribu- 
tions, since Mn ~ N c and L s ~ V-^c> ^ns ~ V-^c) as follows from the analysis of 
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twist -2 parton distributions in the large- N c limit in Ref.[Q- Thus, it makes sense to 
consider both kinds of corrections simultaneously in the large-iV c limit. 

We now compute the isoscalar part, B$, which is leading in the l/iV c -expansion. 
Since Mp oc (p/R) 2 , it is sufficient to consider only the leading contribution in Mp - 
the leading "ultraviolet divergence", which is contained in the lowest-order term in the 
expansion of the matrix element in gradients of the classical pion field of the nucleon 
(see Ref. |28| for details). To derive the gradient expansion we express the nucleon matrix 
element of Eq.( |3.22D in terms of the Euclidean Green function of quarks in the background 
pion field, cf. Eq.fl3.lp, 



G(x,x 4 ;y,y 4 ) = (x,x 4 | 



^.d-tMFi-d^U^Fi-d 2 ) 



|y,2/4>, 



(3.28) 



where |x, x 4 ) denotes a set of formal "position eigenstates" . After integration over the 
soliton center in order to project on the nucleon state with zero three-momentum, and 
over soliton orientations in isospin space, one has (Euclidean conventions) 



Be 



N r 



M 



N 



n a np 



d x tr 



- ll . d - lM F{-d 2 )U^ sss F{-d 2 ) 



(x, 0| F'\-d 2 ) B(-d 2 ) d 2 id^p 
x,0) 



(3.29) 



where n a are the Euclidean components of the light-like vector, and the trace runs over 
Dirac spinor and flavor indices. Expanding in gradients of the classical pion field, and 
evaluating the R.H.S. of Eg. ( pT29|) using plane-wave states, we find in leading order (i, j = 
1,...3): 



Be 



2N c M 2 h 

M N 
2N C M 2 I 2 



riiUj 



3M 



N 



J d z x tr[$l&Jx)0 i l7 dl-1I (x)] 

d 3 x tr[d l U i class (x)d l U class (x)}. 



(3.30) 



In the last step we have averaged over the arbitrary orientation of the three-dimensional 
vector, n: n,:n, — > (%/3. Here, I 2 denotes a quark loop integral of dimension (mass) 2 : 



l/ t IVj 



x 



d 4 p B(p 2 )F(p 2 ) 
(2tt) 4 f 

F 2 (p 2 ) - 2p 2 F(p 2 )F'(p 2 



3 P 



V2 



F 2 (p 2 ) + F(p)F (p 2 



(3.31) 



where p is Euclidean (p 2 > 0) and F' = dF/d(p 2 ). This integral is a measure of the 
"average virtuality" of the quarks and antiquarks in the nucleon, as determined by the 
instanton-induced form factors, F(p). Parametrically, I 2 ~ p~ 2 , i.e., the average virtual- 
ity is of the order of the square of the ultraviolet cutoff, p~ x . Note that in the expression 
Eq.(|3.3ip we have set M = and kept only the leading contribution in Mp. 
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The result for the gradient expansion of Bs, Eq.( |3.30| ), should be compared to that 
the matrix element of the twist-2 spin-2 operator, Eq. fl2.9|) : 



L p) = ^!^| d 3 xtr[5i[/ t ia£s(x)5i[/dass(x)]) (332) 



IN 

where Jo is a loop integral, which is parametrically of order log(Mp): 

f d*p F 2 (p 2 ){F 2 (p 2 )-2p 2 F(p 2 )^(p 2 ) + |p 4 [F ,2 (p 2 ) + F(p)F' , (p 2 ) jj 
~ J (2tt) 4 [p 2 + M 2 F 4 (p 2 )] 2 " ( ' 

Up to terms proportional to derivatives of the form factors it can be identified with the 
pion decay constant in the effective chiral theory |TT|: 



rvj 



AM 2 N C I (3.34) 

[F% = 93MeV). Comparing Eq. (|3~30|) with Eq. (|3~32|) we obtain 

S '' h ^ MH \ (3.35) 



if h f;- 



Numerically we find I2 = 0.0053 p . Using the standard parameters for the average 



instanton size and the dynamical quark mass, p _1 = 600 MeV and M = 350 MeV (TT 
and setting L^ equal to unity, which is the result obtained in the instanton vacuum in 
leading order of p/R f28[, we obtain from Eq.( |3.35|) a value 



B s = 0.9 p' 1 = (570 MeV) 2 . (3.36) 

Thus, Bs has a large numerical value due to the smallness of the instanton size. In 
particular, Bs is of order unity in the packing fraction of the instanton medium, 

Bs ~ p- 2 ~ (|)°- (3.37) 

A comment is in order here concerning the accuracy of the gradient expansion of the 
nucleon matrix elements, which takes into account the leading contributions in Mp. It is 

known that in the case of the twist-2 spin-2 operator the gradient expansion, Eq. (|3.32|) , 

(2) 
underestimates the result of the exact calculation of the matrix element, L s = 1, by a 

factor of 2/3 (see Ref.[|£| for a detailed discussion). This is apparent when one compares 



Eq.( |3.32|) with the well-known expression for the nucleon mass in gradient expansion 

M N = ^ |rf 3 xtr[^f/ c t lass (x)^f/ class (x)]. (3.38) 

The reason for this "2/3-paradox" is that the momentum sum rule, L s = 1, requires 
the equations of motion for the classical pion field of the nucleon to be satisfied; however, 
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the latter have a non-trivial solution only if finite terms in Mp are taken into account. 
An analogous factor of 2/3 occurs in the gradient expansion expression for B$, Eq. (|3.30|) . 



If we had compared our result Eq.( |3.30| ) to the expression for nucleon mass in gradient 



expansion, Eq. (|3.38|) , rather than the twist-2 matrix element, Eq.( |3.32|) , we would have 
obtained an estimate for B$ smaller by a factor of 2/3. We believe the comparison with 
the matrix element of the twist-2 spin-2 operator to give the more realistic numerical 
estimate. Note that this question can ultimately be decided only in a theory which 
correctly takes into account all orders of Mp, that is, of the instanton packing fraction, 
p/R. Since our effective operator has been derived only in the leading order of p/R, it 
would not be sensible to "improve" the calculation of its nucleon matrix element beyond 
this accuracy. At the present level, the difference between the two estimates should be 
seen as indicating the theoretical uncertainty of our result. 

In addition, we have estimated the nucleon matrix element Bs by computing separately 
the contribution of the bound-state level and the Dirac continuum of quarks to the nucleon 
matrix element. The former can be computed exactly, using the known bound-state 
level wave function; the latter can be estimated using an "interpolation formula" (see 
Ref.||16|), which becomes exact in various limiting cases and correctly reproduces the 
leading ultraviolet divergence given by Eq. fl3.30p . The numerical result lies between the 



two possible values obtained from gradient expansion, so all estimates are well consistent. 

3.4 Nucleon matrix elements of the four— fermionic operators 

We now compute the nucleon matrix element of the four-fermionic ( "diquark" ) twist-4 
operator, C a p, Eq.( |2.22j ), which arises from contributions to the DIS cross section of the 
type shown in Fig.|I] (b). 

First we note that the four-fermionic twist-4 operator, C a /3, Eq. fl2.22|) , being a local 
operator, does not have contractions with closed loops which give rise to parametrically 
large (~ p~ 2 ) integrals; the corresponding loop integrals are of type / d 4 k (k^/k 2 ) and 
vanish because of rotational invariance. This is in contrast to the effective four-fermionic 



operator obtained from the quark-gluon operator, B a p, Eq.( |3.12| ), which has a "quadrat- 



ically divergent" contraction, see Fig|J (c). Consequently, the four-fermionic twist-4 
operator, C a p, cannot be represented by a one-body quark operator of the type Eq.( |3.17| ) 
in the effective theory. Its matrix element in a hadron is not related to the virtuality of 
quarks in the bound state; it rather measures two-body correlations between quarks. For 
this reason we expect the nucleon matrix element of C a p to be not of order p~ 2 , but of 
order M 2 , i.e., parametrically small. 

We are interested in the isosinglet matrix element of the operator Eq. (|2.22|) , i.e., the 
average of proton and neutron matrix elements. Because of isospin invariance 

^E e / e /' [(P\ c *p,ff>\v) + (n|C a/ a,//'|n)] 

x x 

= j(e u + e d ) 2 (p\C Qf3:S s\p) + jtt (e u -e d ) 2 (p\C af3yv \p), (3.39) 
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where 

(\a \ f - X a \ 

Z^/yTaTsV'/l IE^/'Y^75^/' I , (3.40) 

H^/ r6 y7a75^/J (5^/'^y7/375^/'J ( 3 - 41 ) 

are manifestly flavor-singlet operators (the sum over a = 1, ... 8 and b = 1, ... 3 is im- 
plied). 

Since g 2 ~ 1/-/V C one finds that the reduced nucleon matrix elements [cf. Eq. (|2.18|) 1 



of the four-fermionic operator, C$s an d CW> are of order N®, and thus of the same 
order as the matrix elements of the twist-4 quark-gluon operators. The nucleon matrix 
elements can be computed in the chiral soliton picture of the nucleon using standard 
techniques. After integrating over soliton translations and rotations, we again express the 
nucleon matrix elements in terms of the quark Green function in the background pion 
field, Eq.( p.28| ) (Euclidean conventions): 

°ss 1 g 2 N 2 f * . ™ n n J * 1 „ n J 1 



n a n p J d 3 xtr[G(x,0;x,0) < fe > 7 Q 7 5 G(x, 0; x, 0) < ft > 7/975]- 

(3.42) 



C vv f 2M N j | r- 



Expanding in derivatives of the pion field we now find in leading order 

C S s \ 4g 2 N 2 M 4 ( r d 4 k F\k) \(t d A k' F\k') 



C vv \ 3Mjv \J (2Tty[k 2 + M 2 F\k)} 2 ) \J (2ixy[k' 2 + M 2 F\k') 

+ 



2 



X 



J d 3 x tr[5 i C/ c t lass (x)9 i C/ class (x)]. (3.43) 



The Euclidean momentum integrals here are parametrically of order log(Mp) ("logarith- 
mically divergent") and, up to terms involving derivatives of the form factors F(k), which 
are not essential here, can be related to the pion decay constant, see Subsection p73| . Com- 
paring Eq. (|3.43|) with the gradient expansion result for the matrix element of the twist— 2 
spin-2 operator, L s , Eq. (|3.32|) , we find 



^,2 z?2 

r - r - 9 tJL t^ 



(3.44) 



Thus, parametrically 



Css. (\x - :i/ J i(»R(3//r) - (|) . (3.45) 
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(2) 



Inserting Eq.( |3.44| ) in Eq.( |3.39| ), taking into account the quark charges, and using Lg = 1, 



we finally obtain 

^Ee /e/ ,[(p|C aA// ,|p) + (n|C aA// ,|n)] = ^A « (50MeV) 2 . (3.46) 

In the last step we have taken a value of g 2 at the scale p~ l = 600 MeV of g 2 ~ 8.8, 
corresponding to /3 = 8n 2 /g 2 = 9.0 [15|. Thus the matrix element of the four-fermionic 



operator, C a/ 3, is significantly smaller than that of the quark-gluon operator, B a p, in spite 
of the fact that the coupling constant at the low scale is rather sizable. The reason for 
this is the different parametric order of the matrix elements: M 2 for the four-fermionic 
operator vs. p~ 2 for the quark-gluon operator. 



In Subsection [3.2| we saw that in the instanton vacuum the matrix element of the 
quark-gluon operator A a p, Eq. fl2.16j ), is suppressed, since the function of the gauge field 
appearing there is zero in the field of one instanton, cf. Eq.( |3.7| ). It is instructive to 
compute the nucleon matrix element of the four-fermionic operator, Eq. fl2.20 ), which is 
obtained from A a p by applying the QCD equations of motion; this allows us to check the 
consistency of our approximations. The calculation is analogous to that for the operator 
C a p, with only the singlet times singlet structure, Eq.( |3.40| ), present. We find 

^four-fermion = ^4^ ~ M 2 k)g(Mp) ~ (£)\ (3.47) 

which is consistent with the fact that the corresponding quark-gluon operator is zero at 
one-instanton level, Eq. (|3.7|) . With g 2 = 8.8 (see above) and Lg = 1 Eq.( |3.47|) gives 
As, four-fcrmion = (195 MeV) 2 , which is an order of magnitude smaller than the correspond- 
ing matrix element of the "true" quark-gluon operator, B$, Eq. (|3.36|) . This shows that 
the parametric suppression of As relative to Bs is indeed borne out by the numerical 
values. 

4 Comparison with DIS data 

We now want to confront our results for the twist-4 matrix elements with the data on Fl 
and F2 from deep-inelastic scattering at moderate Q 2 . We take the usual "phenomeno- 
logical" attitude towards power corrections, namely, assume that there exists a range of 
Q 2 (typically 1 ... 10 GeV 2 ) sufficiently large so that a perturbative QCD treatment with 
individual power corrections is justified, but sufficiently low for the power corrections to 
be noticeable. We shall not be concerned here with theoretical questions related to the 
principal accuracy of the perturbation series, renormalon ambiguities etc.; for this we 
refer to the literature [31, |32], |33| . 

As was shown in Section RL the instanton vacuum implies that among the twist-4 ma- 
trix elements appearing in the l/Q 2 -corrections to Fl and F 2 , Eqs.( [2.17| ) and ( 2.22J) , the 



one of the quark-gluon operator B a p is parametrically leading and has a large positive 
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value, while those of the operator A a p and the four-fermionic operator C a p are para- 
metrically suppressed. We notice that the dominant matrix element, B, enters with a 
large (and positive) coefficient in the twist-4 corrections to F L , Eq. fl2.15p , but with a 
much smaller coefficient in F 2 , Eq. (|2.21|) . Furthermore, in the case of F 2 , Eq.( |2.21| ), the 



parametrically suppressed matrix elements A and C have comparatively large coefficients. 
In the case of Fl the situation is more fortunate. Not only does the dominant matrix 
element in the instanton vacuum, B, enter with a coefficient of order unity, but also the 
coefficient of the parametrically suppressed matrix element A is small. Thus, for Fl we 
can make a quantitative prediction for the sign and magnitude of the twist-4 contribu- 
tion. In addition, for this structure function the twist-2 contribution is non-zero only in 
order as, cf. Eq. fl2.7j ), so the relative importance of power corrections is larger than in 
F 2 , which is indeed seen in the data for R ||, [|. We shall therefore concentrate on Fl in 
the following. 

We would like to compare our theoretical result for the scale dependence of F L 
with measurements of the ratio R for a deuteron target [|], Q. We assume here that 
^l ~ i( l + ^l)> that the deuteron structure function approximately corresponds to 
the isoscalar combination of proton and neutron structure functions. The full theoretical 
result for the Q 2 -dependence of the second moment of F L consists of the twist-2, target 
mass, and dynamical twist-4 contributions, see Section |2|. We evaluate the twist-2 and 
target mass contributions, Eqs. (|2~7"l ) and (|2.13|) , for the isoscalar target, using the GRV94 



NLO parametrizationQ see Ref . [|34"f for details such as the number of quark flavors, Aqcd, 
etc. For the dynamical twist-4 contribution we take the values obtained from the in- 
stanton vacuum, Eq.( |3.36| ), and neglect the contribution of strange quarks in the twist-4 
corrections. This results in a total twist-4 contribution of 

Mf(2,Q 2 )— « I^ X ^ + e ^ = (26 °g 2 eV)2 - (4- 1 ) 

It is instructive to first compare the relative magnitude of the power corrections and the 
0(as)-perturbative contributions to Ml, see Fig.||. The dashed line shows the dynamical 
twist-4 contributions obtained from the instanton vacuum, Eq.( |4.1| ). The dot-dashed line 
is the sum of twist-4 and target mass corrections, Eq.( [2.13| ), i.e., the total l/Q 2 power 
corrections. The solid line, finally, is the total result, in which also the twist-2 contribu- 
tion has been added; the variation of the twist-2 contribution with Q 2 over the interval 
shown is much weaker than that of the power corrections. One sees that the total power 
correction is sizable, which is clearly a consequence of the missing tree-level [^(a^)-] 
perturbative contribution, and that power corrections are dominated by the dynamical 
twist-4 contribution, which is much larger than the target mass corrections. These facts 
combined make the longitudinal structure function an excellent tool for extracting quan- 
titative information about higher-twist matrix elements. 

We also compare our total theoretical result for Mj* to order l/Q 2 with the moment 
extracted from the fit to the SLAC data for R d [HI, which describes well also the more 



4 The use of a NLO parametrization for the twist-2 parton distributions is consistent with the use of 
the O (as )-expressions for the Wilson coefficients in the twist-2 contribution. 
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recent deuteron data from the NMC experiment |35|| . In order to obtain an estimate for 
isosinglet F d to accuracy 1/Q 2 we have extracted from the fit for R d in Ref.0 the total 
contribution up to order 1/Q 2 , and combined this with the fit to the NMC data for i 7 ^ 
36|, which also includes l/Q 2 -corrections. Substituting the fits for R d and F d in Eq.( [2.4| ), 
retaining only terms up to order 1/Q 2 , we obtain a rough estimate of Fl to order 1/Q 2 . 
To be able to compute the moment we have estimated the contribution of the unmeasured 
small-x region using the GRV parametrization; this contribution turns out to be at the 
few percent level. The "experimental" result for the moment thus obtained is shown by 
the dotted line in Fig|| it should be compared to the solid line giving the total theoretical 
result. As can be seen the agreement between the total results is rather satisfactory. This 
indicates that the twist-4 contribution derived from instantons is of the right order of 
magnitude. 

We note that our results for the twist-4 matrix elements obtained from the instanton 
vacuum are consistent with those of the analysis of power corrections to the second mo- 



ments of Fl and F2 of Ref. [24], which was based on NMC and SLAC data. Assuming only 



that the matrix elements of the four-fermionic operators, A and C, Eq. Q2.16 ) and (|2.22|) 



are not strongly different in magnitude, the authors of Ref. [24"] concluded that it is likely 



that the matrix element E>s has a large positive value of the order of (380 ± 95MeV) 2 
at jj? = 5 GeV 2 , somewhat smaller than our estimatef] We have repeated the analysis 



of Ref.f24] using the results of a recent re-analysis of SLAC and NMC data of Ref. [ 



where the twist-4 corrections were parametrized using a renormalon fit, and found only 



small changes to the values quoted in Ref. [24]. In particular, as explained in Subsection 



the instanton vacuum provides a parametric explanation why the contributions of the 
four-fermionic operators are suppressed. 

A more conclusive comparison with phenomenological power corrections could be made 
if we had results for the x-dependence of the twist-4 contribution. In the OPE approach 
of Refs. m, [U this would require the computation of the nucleon matrix elements of higher- 
spin operators describing power corrections to higher moments, a problem which is cur- 
rently under investigation. In contrast, in the renormalon approach to power corrections 
[2~T| , B7| , |3"8] | the x-dependence of the l/Q 2 -correction is obtained with relative ease, while 
it is impossible to predict the absolute normalization. Fits to the experimental data for 
R and F2 indicate that the renormalon formulas describe the x-dependence of the data 
rather well, at least for large x |6|, |7| . In this situation it is interesting to make a rough es- 
timate of the full x-dependent twist-4 contribution simply combining the two approaches, 
i.e., taking the x-dependence from the renormalon formula and determining the coeffi- 
cient by the instanton vacuum calculation. We use the renormalon expression for the 



twist-4 correction to Fl given in Refs. ]F2T] , pq] , including only the non-singlet coefficients, 
which we evaluate using the GRV94 NLO parametrization for the twist-2 parton distri- 
butions. (The logarithmic Q 2 -dependence of the twist-2 distribution may be neglected 
in this context.) Combining this with the twist-2 plus renormalon fit to F 2 of Ref.|J 
we obtain an estimate for R(x, Q 2 ) to order 1/Q 2 , shown by the solid line in Fig.|j. It is 

5 Here we have translated the results of Ref. [E4| into our conventions. 
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Source 


K 2 /GeV 2 


Scale / GeV 2 


Present work (instanton calculation) 

Choe et al. (data analysis) |2~3|] 

Yang and Bodek (renormalon fit) || 

Ricco et al. (renormalon fit) |7J 

Stein et al. (MS renormalon model) (21], 22 


0.068 

0.030 ±0.015 
0.034 ±0.002 
0.060 ±0.020 
±0.032 


-0.36 

5 
1 
1 
1 



Table 1: The total twist-4 contribution to the second moment of Fl, parametrized here 
as Ml(2, Q 2 ) twlst ~ 4 — k 2 /Q 2 , as extracted from various fits to experimental data or model 
calculations. 



compared to the recent NMC data for R(x, Q 2 ) [|35| , taken at various Q 2 . In view of the 
simplicity of this model the agreement is quite reasonable. Note that this model is not 
applicable in the small-x region. 

Finally, in Table [l| we compare the instanton vacuum result for the total twist-4 
contribution to the second moment of Fl, Eq.( [4.1|) , with results of other approaches. In 
the analysis of Ref.|| the higher-twist contributions to R and F% at large x were modeled 
by the non-singlet renormalon contribution p2| [37], |38|, with the coefficients treated as a 
fitting parameter. Using these fits to construct Fl, and computing the second moment, we 
obtain the value quoted in Table |l|. Also given is the result of a recent similar renormalon 
analysis J7|. Both values are in reasonable agreement with the instanton calculation. 
Finally, within a given renormalization scheme it is possible to determine the absolute 
magnitude (but not the sign) of the renormalon contribution; we give in Table |I] the value 
obtained in Ref.[^2j in the MS scheme. 

In the so-called transverse basis of twist-4 operators the twist-4 contribution to F L 
can be related to the "average transverse momentum" of the quarks in the nucleon ||. 
Forms such as 



Mt(2,Q 2 ) 



2\ twist— 4 



/ 2 2\ \ TV 

2 {6u + ed) ~Q 2 ~ 



(4.2) 



are also frequently used as a phenomenological parametrization of the twist-4 contribution 
to Fl. Comparing Eq. fl4.2|) with the OPE result in the "longitudinal" basis, Eq. (|2.15|) , 
and keeping only the dominant matrix element in the instanton vacuum, Bs, we see that 
our result for B s , Eq.( |3.36j ), corresponds to 



(4) 



16 



It is interesting that (kj) is parametrically of order p 



:B S « (250 MeV) 2 . (4.3) 

i.e., it is determined by the 
instanton size, or, in other words, the "size" of the constituent quark, not by the size of 
the nucleon, which is of order M _1 . Thus, the smallness of the instantons explains the 
relatively large value of (k 2 -,) found in power corrections. 
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5 Conclusions 

In this paper we have studied the twist-4 matrix elements describing l/Q 2 -corrections 
to unpolarized structure functions in the instanton vacuum. We have shown that the 
dominant twist-4 matrix elements are those of the "true" quark-gluon operators, i.e., 
where the gluon field cannot be eliminated using the equations of motion, not those of 
the four-fermionic (or "diquark" ) operators. The typical scale for the matrix elements of 
the former is the square of the instanton size, p~ 2 , while the latter are of order of the dy- 
namical quark mass squared, M 2 , and thus suppressed by the instanton packing fraction, 
(p/ R) A . This result is contrary to the widely held belief, based on a naive "constituent 
quark" picture, that it is the diquark correlations which are mostly responsible for power 
corrections to structure functions. Rather, in our picture, the dominant power corrections 
come from the interaction of the quarks with non-perturbative vacuum fluctuations of the 
gauge fields of small size. 

The numerical values obtained for the various twist-4 matrix elements follow their 
parametric ordering. The dominant quark-gluon operator, B a p, has a large positive value, 
while the matrix elements of the diquark operators are an order of magnitude smaller. 
This seems to be in qualitative agreement with the results of a phenomenological analysis 
of power corrections to F% and R d . We are not claiming, of course, that the observed 
power corrections could only be explained by instantons. Contrary to the phenomenon of 
chiral symmetry breaking, which can directly be linked to the presence of topologically 
non-trivial vacuum fluctuations, in the case of power corrections to structure functions 
there is no a priori reason why instantons should play a role. Nevertheless, it seems 
that the basic assumptions made in the instanton vacuum, namely the dominance of 
instanton-type fluctuations and diluteness of the instanton medium, are consistent with 
the information we have on power corrections to structure functions. 

We have shown that in the instanton vacuum the effective operator for the leading 
quark-gluon operator B a p can be reduced to a one-body quark operator measuring the 
"average virtuality" of the quarks in the nucleon. The sign and magnitude of the twist-4 
matrix element B depend not only on our picture of non-perturbative vacuum fluctua- 
tions, but also on our description of the nucleon. In this respect it is crucial that the 
picture of the nucleon as a chiral soliton, which is a fully relativistic and field-theoretic 
treatment, is obtained from the effective low-energy theory without making any addi- 
tional approximations. The correct sign of the matrix element B$ is obtained in this 
approach due to the fact that the average momentum of the quarks in the nucleon is 
space-like, which would not be the case, say, in a simple constituent quark model. 

The approach taken in the present work can be extended in many ways. A straightfor- 
ward continuation would be the computation of the matrix elements of twist-4 operators 
giving power corrections to higher moments of the structure functions, and the restoration 
of the ^-dependence of the twist-4 contribution. It is likely that twist-4 QCD operators 
of spin > 2 can again be represented in the effective theory by simple one-body quark 
operators of the form of twist-2 operators with extra derivatives, as in the case of spin 2, 
Eq.(|3.22|). A much more challenging task would be to depart from the notion of power 
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corrections altogether and develop a uniform description of the Q 2 -dependence of struc- 
ture functions, valid also at small Q 2 , e.g. in the spirit of Ref.||19||. The instanton vacuum, 



which on one hand has a clear relation to perturbative QCD, on the other hand describes 
non-perturbative phenomena in the infrared domain, may be a good starting point for 
such attempts. 

Finally, we note that the methods developed here can be applied to study also power 
corrections to exclusive processes such as deeply-virtual Compton scattering and hard 
meson production, which may be of great practical importance |39[]. For these processes, 
however, up to now no systematic classification of power corrections in QCD is available. 
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A Operator product expansion for unpolarized DIS 

In this appendix we present a short derivation of the QCD predictions for the (^-depen- 
dence of the second moments of F 2 and Fl using the operator product expansion. 

From the definition of the hadronic tensor, Eq.( |2.1| ), explicit expressions for the struc- 



ture functions F 2 and Fl can be obtained by contracting Eq.(2.1) with p' l p l/ , taking its 
trace, and solving the resulting equations: 

F L (x,Q 2 ) = ^ (M 2 N W + 2W PP ) + O ( ±) , (A.l) 



Q2 V - iv --wj - y Q4 



F 2 (x,Q 2 ) = -xW + (M 2 N W + W PP ) + 0( — \, (A.2) 



where 



W = W* W pp = pV^- (A.3) 



We have kept only terms up to order 1/Q 2 . Note that W and W pp generally contain terms 
of any order in 1/Q 2 and should be substituted in Eqs.(SH, |A.2| ) consistently. 



The QCD description makes use of the fact that, by the optical theorem, the hadronic 
tensor is related to the imaginary part of the forward virtual Compton amplitude, 

T pu (p,q) = i f d^ze^ z ^\T{J^z)J v {Q)}\p)- ( A -4) 

As a function of the photon energy, q°, it exhibits a cut in the inelastic region. To express 
the hadronic tensor in terms of the Compton amplitude it is convenient to regard the 
latter as a function of the dimensionless variable 

U = x = ~W' (A ' 5) 

in this variable the cut runs along the real axis from —00 to —1 and from +1 to +00. In 
the case of electro-magnetic scattering the Compton amplitude is crossing-even, and thus 
an even function of to. The relation to the hadronic tensor is provided by the formula 

T^ ) , /-i . n _ x f W(x 



4 / dxx n - 1 { , (A.6) 



Tpp ;U) n I Jo I W pp [x 



where T and T pp are defined in analogy to Eq.( |A.3| ), and T^n , T PPA ,n denote the coefficients 



of uf 1 in the power series expansions of T and T pp in u in the domain |u;| < 1. In particular, 
for the second moments of Fl and F 2 one obtains in this way: 

jT 1 dxF L (x, Q 2 ) = ± (M 2 N T^ + 2T pp ^) + O (±A , (A.7) 

J Q 1 dxF 2 (x,Q 2 ) = ~T^ + -L(m 2 n T^ + 3T pp ^) + oUyj. (A.8) 
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The general strategy is now clear. One computes the virtual Compton amplitude in the 
Bjorken limit, Q 2 — > oo and uo fixed, expanding in 1/Q 2 , and pick up the coefficients in T 
and T pp of the given powers in u>. 

The dominant contribution to the virtual Compton amplitude in the Bjorken limit is 
given by the "handbag" diagram describing the scattering of the photon off a single quark. 
At level 1/Q 2 one must consider generally two types of contribution, corresponding to the 
different contractions of the quark fields in Eq.(A.4) shown in Fig.|T]. Contributions of type 
(b) can be computed using standard Feynman diagrams and give rise to matrix elements 
of four-fermionic operators in the nucleon (see below). Contributions of type (a) require 
the calculation of the quark Green function in the background of the non-perturbative 
gluon field of the nucleon. This can conveniently be done using the Schwinger method 
|I|, [H]]. In this approach the contribution (a) to the Compton amplitude is given by the 
formal expression (q = q ar y a etc.) 



T,u = E4W/ 



-7m 



1 



P + q 



lv-lv 



P-q 



In 



^/|P>> 



where 



P M = iV M = idn + Ap, [P^Pu] = iF^. 

Eq. (|A.9|) is understood as a starting point for expanding in local operators: 
1 q qPq qPqPq qPqPqPq 

P~^ q = f'T^f Iff ¥)^ + '" 



(A.9) 



(A.10) 



(A.ll) 



To order 1/Q 2 all contributions to the second moments of Fl and F 2 , Eqs.( |A.7| ) and 
( |A.8| ), are contained in the O(P) and 0(P 3 ) terms of Eq.( |A.li|) . We now list these 
contributions. For simplicity we suppress the sum over flavors and the quark charges; 
they will be inserted again at the end of the calculation. 

The 0(P)-term of Eq.( A.ll ) contributes to T^i and gives rise to the scaling part of 
the second moment of F 2 : 



(Q 



T (order P) = T/ (order P) 

Using the identity 

qP = 2q-P-Pq 
as well as 7^7^ = — 2q, one can write this as 



2^2 



{p\^qPqi^ip\p). 



(A.12) 



(A.13) 



T(order P) 



(q 



2\2 



q a q^\^P a l^\p) + 



(A. 14) 
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where we have dropped terms with contracted g's which do not contribute to the cu 2 -term. 
We see that this contribution is entirely given by the nucleon matrix element of the spin-2 
twist -2 operator, 

(p\$P al ^\p) = 2L { - 2 \p a p (3 -\p 2 g a( ,)- (A.15) 

the traceless structure here is a consequence of the QCD equations of motion: ifjPip = 0. 
Substituting Eq. (|A.15| ) in Eq.( |A.14|) one finds 



T(orderP) = -4L (2) a; 2 + . . . , (A.16) 

and from Eq.( |A.8| ) one obtains 

£dxF 2 (x) = L^+o(-^\. (A.17) 



The 0(P 3 )-term in the expansion Eq.( A.li ) contributes to the u 2 as well as the de- 
part of T, and to the u; 4 -part of T pp . It gives rise to the dynamical twist-4 as well as the 
target mass corrections to the second moments. One can easily see that the dynamical 
twist-4 contributions are contained in Tji and T pPtUJ 4, since a twist -4 matrix element can 
contribute at most two powers of the target momentum. (In particular, this means that 
the twist-4 contribution to the second moment of Fl comes only from T ppA ,4.) Target 
mass corrections, on the other hand, arise from all terms. We begin with 

T pp (orderP 3 ) = -^-pV {p\^^qPqPqP qi ^\p). (A.18) 

Making repeated use of the identity Eq.( [A.13|) , dropping along the way all terms with 
qq = q 2 in the numerator which do not contribute to the u; 4 -piece, one can reduce this to 

T pp (orderP 3 ) = ™pV (p|$(2g ■ P)\qj„ip\p) + • • • 

16 

(2g-p)gVgV-M^VgV (plW^iVyWb) + ••■ (A.19) 



2^4 



The reduction of the trace term, T, is slightly more complicated, since one has to retain 
the g 2 -terms, which can contribute to the u 2 -piece. Again using Eq.( |A.13| ) one first 
obtains 

2 - - ~ * 

T(order P 3 ) = ■p^j{p\ip r y fJ 'qPqPqPq r Y fi ip\p) 

= -^<p| {(2q ■ P) V^ ~ Q 2 [(2q ■ P)-fPl^q • P) 

+ {2q ■ P) 7 ^PP 7m + ^PP qill (2q ■ P)] } \p) + . . . (A.20) 
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Anticommuting now the matrices 7 M from the left to the right position, and taking note 
of the QCD equations of motion 

ip...Pif> = tpP...ip = 0, (A.21) 

one finally obtains 

Tforder P 3 ) = -^- \-2q a q p q<q & + q 2 (-q a q^g pS + q a q S g Pl + q 1 q & g aP ] 
(g J ) 4 L v yj 

x{p\iJ;P a P P 7l5 i;\p) + ... (A.22) 

We see that both Eqs.( |A.19| ) and ( |A.22|) are determined by the matrix element of the 
reducible rank-4 tensor operator, ipP a PpP 1 '-fsifj. Only the part symmetric in the indices 
a and 7 enters. On general grounds, this matrix element can be parameterized asQ 

(p\^P {a} PpP {j} ^\p) = 2L®S alh6 

+ 2&! (g a pS^ + g^pSas) + 2b 2 g ai Sps 

+ 26 3 {ga&Sp-f + g-fsSap) + Zb/tgpsSa-f 

+ terms g a pg^s etc. (A. 23) 

Here S a p__, denotes the totally symmetric traceless tensor constructed from the vector p^. 

Sa/3^5 = PaPpP-yPS ~ \p 2 (gapPyPS + • • • + g-ySPaPp) 

+i$P 4 {9a09fS + da^gps + gasgp-y) , (A.24) 

S a p = PaPp ~ \p 2 g a p- (A.25) 

The coefficients in the ansatz Eq.( |A.23|) can be expressed in terms of matrix elements of 
operators of definite spin (and twist). L^ is determined by the matrix element of the 
unique spin-4, twist-2 operator 

{p\$P {a PpP ll5} iP\p)- traces = 2L^S a p l5 . (A.26) 

When inserted in Eqs. ( |A.19| ) and (A.22) this part of the matrix element gives contributions 

T(orderP 3 ) = L^ (-4u 4 - ^^u 2 ) , (A.27) 

T pp (orderP 3 ) = ^L^w 4 , (A.28) 

and thus gives rise to the "target mass" corrections!] 

r 1 T^ Mir 

/ dxP 2 (x,Q 2 ) targctmass = -r--^, (A.29) 

Jo 2 Q z 



/ dxP L (x,Q 2 ) targctmass = -3L^-^-. (A.30) 

Jo Q 2 



6 On grounds of Lorentz invariance and symmetry in a, 7 Eq.( A.25| ) could contain also a structure 
i{^<ri3-ysS aa + to-pasSa-f), which, however, is not related to the matrix element of a twist-4 operator. 
7 This result agrees with Ref. [pi when one expands the Nachtmann moments there in Mf^/Q 2 . 
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The coefficients h\ ... 64 in Eq. ( |A.23| ) can be expressed in terms of the matrix elements 
of certain spin-2 twist-4 operators related to contractions of the rank-4 operator. We 
consider the set of operators 



A aP = ^ a [V\ F 7/3 ] ^ = ^ a {D^F^, 



(A.31) 
(A.32) 
(A.33) 



where symmetrization in a, (3 and subtraction of traces is implied. Their matrix elements 
are given by 



(p\Aa(}\p) 



2AS a , 



etc. 



(A.34) 



The expressions for b\ ... 64 can be found by taking traces of both sides of Eqs. (|A.23|) . By 
straightforward calculation one finds 



a/3 



7 a 7 



,aS 



> (p\ijP {a} PpP {7} -fs^\p) 



g 



as 



{ (D-B) S 7S , 
(D-B + A)Sp S , 
0, 
D Sap. 



(A.35) 



The operator B a p comes into play when writing P^P^ = PP — (i/4)[7 P , 7 CT ]F pCT and 
making use of the three-gamma identity 



7a7/37 7 = fi'a/377 - 9ajlP + 9fala ~ ICa^Sl 75- 

We are using the conventions 



75 



-zWtV* 



tr 



VWtV 



IE 



a/375 



,0123 



(A.36) 



1. (A.37) 



Solving the system of equations obtained by contracting Eqs. ( |A.23j ) for 61 ... 64 one ob- 
tainsQ 



h 

b 2 
h 

by 



1 
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-A -3B +2L>, 

5A SB +4D, 

-A +B -2D, 

A +B +AD. 



(A.38) 



8 This result is in agreement with that of Ref.|l|. Note that in Ref.jlj t a = X a , and that in Ref. [Q the 
operator A^ v is defined with different sign and relative coefficient. 
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Inserting the spin-2, twist-4 part of the matrix element in Eqs.( [A.19|) and ( |A.22| ) one 
obtains 

T(orderP 3 ) = — (A + 5B)cu 2 + . . . , (A.39) 

T pp (orderP 3 ) = -(-A + 3B)u* + . . . , (A.40) 



and thus from Eqs.( |A.7[ ) and ( |A.8| 



'dxF^Q 2 )^- 4 = -L(-A + 3B), (A.41) 

o 4Q Z 

f 1 dxF 2 (x,Q 2 ) t - ist - i = -L(-5A-B). (A.42) 

As noted in Ref . |l| the operator D pv does not appear in the final result for the dynamical 
twist-4 contributions. 

To obtain the final expressions one still has to take into account the quark charges, 
i.e., replace in all the above matrix elements of quark bilinear operators 

^...v -► J2 e }$f---^f- ( A - 43 ) 

B Minkowskian vs. Euclidean theory 

Here we state our conventions for passing from the Minkowskian to the Euclidean field 
theory, in which the instanton vacuum is defined. After continuing the time dependence 
of the fields to imaginary times we introduce Euclidean vector components according to 

(x°)m = -i(xi)E, (x 1 )m = (Xi)E- (B.l) 

The corresponding rules for derivatives are 

(<%)m = «(<9 4 )e, (di) M = (di) E , (B.2) 

and the Euclidean gauge potential is introduced analogously, 

(A))m = KA^e, (Ai) M = (Ai) E , (B.3) 

so that the relative sign between derivative and gauge field in the covariant derivative, 
Eq.(|2.10D, remains the same in the Euclidean theory. The Euclidean field strength is 



again defined by Eq. (|2.11|) . The dual field strength in the Euclidean theory is defined as 
(F^) E = (l/2)e flupa (F pa ) E , with e 123 4 = 1. 

For use in the Euclidean theory we introduce hermitean gamma matrices according to 

(7°)m = (74)e, (Y)m = z(7i) E , (B.4) 

and use if>* = iip for the Dirac conjugate fermion field (not to be confused with the 
hermitean conjugate of ip). 
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(a) 



(b) 



Figure 1: The two types of twist-4 contributions to the DIS cross section, (a) Interference 
between scattering of the photon off a free quark and scattering off a quark interacting 
with the non-perturbative gluon field of the nucleon. (b) Interference between scattering 
off two different quarks in the nucleon, with a perturbative gluon exchange. 
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(a) 



(b) 





(c) 



Figure 2: Graphical illustration of the effective operator approach, (a) The QCD quark- 
gluon operator B a p, Eq.(|3.9|), normalized at // = p _1 = 600 MeV. (b) The effective four- 
fermion operator, Eq. (|3.12| ). The gluon field has been replaced by the field of an (anti-) 
instanton, which couples to quarks through the zero modes, (c) The two contractions of 
the four-fermionic operator giving rise to parametrically large contributions. The loop 



integrals here are "quadratically divergent" (~ p 
defines the effective one-body operator, Eq.( |3.f3| ). 



^-2 



). The sum of the two contractions 
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Figure 3: The (3 2 -dependence of the second moment of Fl for the deuteron, M||(2, Q 2 ), 
Eq. (|2.6|) . Dashed line: Dynamical twist-4 contribution obtained from the instanton vac- 
uum, Eq.( [4.1[ ). Dot-dashed line: Twist-4 contribution plus target mass corrections, 
Eq. ( |2.13D . Solid line: Total theoretical result, including also the twist-2 radiative correc- 
tions, Eq.( p.7|) . Squares: "Experimental" result for Ml(2,Q 2 ) obtained by combining the 
fit of F 2 d of Ref.|36| with the fit of R d of Ref. Q, keeping only terms up to order 1/Q 2 (see 
the text for details). 
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Figure 4: The ratio R d = crf/a^ obtained from the renormalon model (solid line), com- 
pared with the NMC data from Ref.|3S|. Here the x-dependence of the l/Q 2 -contribution 
is predicted by the renormalon formula of Refs. |37|, [38| , while the coefficient is determined 
from the instanton vacuum result, Eq.( (4.1| ). In the data set the values of Q 2 vary with 
x; see Ref. |35| for details. Note that the simple non-singlet renormalon formula is not 
applicable at small x. 
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